Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1359 1768 (N39W51 L=199)  1557    295        4 (1222/1707-1472/1815)  -3.5 W35-W68 N37-N40
      2  660 1408 (N34E15 L=133)   434    108        1 ( 610/1389- 728/1427) -18.4 E18-E09 N31-N35
      3 1756 1332 (N03W56 L=204)   432     63        0 (1721/1311-1786/1350) -31.0 W53-W59 N01-N04
      4  887 1309 (N22E02 L=146)   237     77        0 ( 847/1293- 923/1323)   0.0 E04-W00 N21-N23
      5 1372 1253 (N08W25 L=173)   299     63        0 (1350/1220-1395/1286) -78.7 W24-W26 N05-N10
      6 1810 1172 (S07W56 L=204)   478     27        0 (1787/1152-1835/1197)   0.0 W53-W58 S08-S05
      7  117 1149 (N29E73 L= 74)   741     91        0 (  96/1110- 141/1192) -21.0 E82-E67 N26-N31
      8  934 1165 (N13E02 L=145)   225    120        0 ( 877/1148- 991/1185) -18.4 E05-W01 N12-N14
      9 1857 1077 (S14W59 L=207)  2511    206        9 (1802/ 997-1908/1171)  29.7 W51-W68 S17-S09
     10 1502 1128 (S02W30 L=178)   367    170        1 (1432/1097-1559/1169)  -5.7 W25-W35 S03-N00
     11  390 1051 (N18E40 L=108)  1497    347        5 ( 243/1006- 514/1142)   7.1 E53-E29 N15-N21
     12 1616 1028 (S10W36 L=184)   368     83        1 (1602/ 985-1632/1079)  68.2 W34-W37 S12-S07
     13 1717  882 (S21W42 L=189)   502    110        0 (1696/ 832-1733/ 943)  71.6 W39-W43 S24-S17
     14 1565  796 (S23W29 L=176)   918    202        1 (1495/ 714-1616/ 883)  36.0 W22-W33 S26-S18
     15 1377  742 (S21W15 L=162)   807    278        0 (1332/ 627-1437/ 865)  47.5 W09-W20 S27-S15
     16 1598  688 (S30W30 L=178)   456    120        0 (1546/ 646-1640/ 740)  18.4 W25-W34 S32-S27
     17 1007  536 (S25E13 L=135)   287     87        0 ( 977/ 495-1037/ 573)  38.7 E15-E09 S27-S23
     18  706  351 (S30E41 L=107)   905    253        2 ( 685/ 232- 730/ 474)  43.4 E51-E33 S39-S22
  

    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