Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1512 1585 (N23W51 L= 34)  1128    218        1 (1483/1477-1537/1664)  37.6 W42-W57 N17-N27
      2 1271 1490 (N25W30 L= 13)   384    133        0 (1207/1469-1328/1514) -45.0 W26-W32 N22-N28
      3  220 1391 (N49E67 L=276)   594     91        1 ( 194/1347- 255/1428) -71.6 E72-E59 N45-N51
      4 1153 1378 (N23W18 L=  1)   774    279        1 (1029/1328-1253/1418) -28.3 W09-W22 N17-N26
      5  939 1305 (N25W03 L=346)   235    107        0 ( 917/1269- 968/1354)  31.0 W00-W06 N22-N26
      6  935 1192 (N18E01 L=342)   835    335        3 ( 850/1063-1054/1306) -77.5 E04-W02 N08-N26
      7  384 1248 (N36E39 L=304)   998    177        3 ( 316/1201- 464/1293) -25.0 E46-E31 N32-N39
      8  533 1135 (N25E28 L=315)   307     76        1 ( 504/1110- 563/1173) -63.4 E29-E25 N23-N27
      9  159 1016 (N25E67 L=276)   778     74        2 ( 142/ 979- 174/1059)  76.0 E70-E64 N23-N27
     10  187  931 (N20E64 L=279)  2358    194        8 ( 133/ 848- 231/1035)  48.0 E72-E57 N14-N25
     11 1429  869 (S13W19 L=  2)   974    313        0 (1365/ 742-1488/1008)  42.7 W12-W25 S18-S06
     12 1714  899 (S21W40 L= 23)   750    167        2 (1688/ 832-1734/ 990)  60.9 W37-W43 S24-S15
     13  668  897 (N08E24 L=319)   572    174        2 ( 608/ 835- 727/ 964) -68.2 E26-E21 N03-N13
     14  373  879 (N14E45 L=298)   269     77        0 ( 326/ 868- 402/ 895)  -8.1 E49-E42 N12-N14
     15 1426  597 (S30W13 L=356)   697    176        1 (1373/ 532-1475/ 665)  31.0 W07-W17 S32-S26
     16 1179  616 (S21E02 L=341)   372    139        1 (1144/ 590-1236/ 661) -45.0 E04-W01 S23-S17
     17 1120  297 (S42E18 L=325)   546    119        1 (1080/ 256-1159/ 349)  26.6 E23-E13 S44-S38
     18  931  257 (S39E34 L=309)   632    121        1 ( 904/ 201- 952/ 317)  51.3 E38-E30 S43-S33
  

    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