Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1667 (N48W01 L=220)   327     47        0 ( 748/1649- 811/1683) -36.9 E02-W02 N46-N49
      2 1606 1627 (N23W59 L=278)  1524    160        7 (1541/1567-1661/1671) -49.4 W54-W67 N19-N26
      3  594 1627 (N51E18 L=202)   395     60        0 ( 562/1601- 621/1661)  26.6 E21-E15 N49-N53
      4 1256 1602 (N31W31 L=250)  2064    419        6 (1064/1544-1428/1643) -28.6 W17-W41 N23-N35
      5 1683 1534 (N16W59 L=279)   992    200        2 (1622/1470-1744/1590)   5.7 W50-W70 N13-N18
      6 1571 1526 (N19W49 L=268)   496    103        2 (1518/1499-1609/1551) -26.6 W45-W52 N16-N20
      7 1659 1339 (N06W47 L=267)  1202    274        4 (1532/1285-1734/1403) -39.8 W40-W53 N02-N12
      8  603 1169 (N20E22 L=198)  1381    398        2 ( 413/1108- 757/1232) -32.5 E35-E13 N13-N27
      9  131 1107 (N25E68 L=152)   719     60        0 ( 101/1080- 166/1135)  -9.5 E73-E61 N23-N27
     10  321  950 (N13E47 L=173)   389     98        2 ( 295/ 918- 357/ 996) -78.7 E48-E43 N10-N15
     11 1463  935 (S11W23 L=243)   277     87        2 (1440/ 910-1499/ 976) -68.2 W22-W25 S13-S08
     12 1399  824 (S16W17 L=237)   317    124        1 (1359/ 789-1451/ 850)   8.1 W14-W21 S18-S15
     13  876  718 (S11E15 L=205)   462    186        0 ( 849/ 662- 948/ 769)  23.2 E17-E10 S14-S08
     14 1875  740 (S35W61 L=280)  1143    135        0 (1816/ 706-1937/ 760)  -4.8 W52-W77 S36-S33
     15 1482  630 (S30W20 L=240)   845    222        0 (1390/ 578-1542/ 709)  18.4 W13-W25 S32-S26
     16  854  439 (S28E25 L=194)  1781    457        2 ( 674/ 361-1055/ 550)   5.4 E41-E09 S29-S24
     17  404  406 (S20E62 L=158)   851     76        4 ( 383/ 362- 423/ 449)  36.9 E67-E58 S23-S17
     18 1230  340 (S44E04 L=215)  1719    359        1 (1062/ 290-1368/ 438)   4.4 E18-W08 S47-S40
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y