Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  958 1540 (N35W09 L=131)   859    240        0 ( 859/1509-1070/1564) -11.3 W01-W16 N31-N36
      2 1590 1486 (N16W49 L=172)   410     46        1 (1566/1465-1619/1506)   0.0 W46-W51 N14-N17
      3 1282 1319 (N14W22 L=144)   358    128        0 (1222/1302-1345/1339) -14.0 W18-W26 N12-N15
      4  284 1259 (N32E49 L= 73)  1880    271        2 ( 198/1160- 381/1333) -32.5 E61-E39 N24-N38
      5  729 1183 (N18E14 L=108)  1502    422        1 ( 549/1122- 898/1292)   2.3 E26-E01 N15-N20
      6  137 1156 (N28E70 L= 53)   837    107        1 ( 114/1116- 156/1196)  76.0 E73-E65 N25-N30
      7 1656  793 (S25W36 L=159)  1266    282        3 (1599/ 675-1703/ 926)  52.4 W30-W41 S31-S18
      8 1373  729 (S22W15 L=137)   614    249        1 (1317/ 619-1421/ 833)  47.7 W08-W19 S27-S16
      9 1856  787 (S31W59 L=181)   900     90        1 (1821/ 752-1896/ 828)  16.7 W53-W64 S32-S29
     10  507  553 (S14E45 L= 78)   672    130        2 ( 487/ 479- 523/ 619)  56.3 E48-E42 S18-S09
     11 1473  583 (S34W20 L=142)   305     41        0 (1443/ 564-1499/ 602)  14.0 W17-W21 S35-S32
     12 1011  512 (S27E13 L=109)   479    103        1 ( 987/ 455-1039/ 558)  50.2 E15-E10 S30-S24
  

    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