Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1252 1814 (N45W54 L= 75)   892    138        2 (1175/1792-1317/1832) -21.0 W47-W62 N42-N48
      2 1280 1457 (N23W29 L= 49)   384     68        2 (1248/1432-1302/1486)   9.5 W25-W31 N21-N25
      3 1088 1426 (N27W16 L= 36)   244     62        0 (1071/1397-1112/1452)  36.9 W14-W18 N25-N28
      4 1795 1319 (S00W61 L= 81)   521     59        0 (1766/1304-1833/1334)  -6.3 W57-W66 S01-N00
      5 1021 1218 (N18W05 L= 25)   744    226        2 ( 972/1128-1091/1301)  28.6 W00-W11 N13-N20
      6  818 1238 (N24E06 L= 14)   213     80        0 ( 781/1222- 864/1258) -31.0 E08-E03 N22-N25
      7  659 1137 (N22E19 L=  2)   586    203        1 ( 563/1087- 733/1168) -35.0 E25-E15 N17-N25
      8  414 1038 (N22E39 L=342)   601    152        4 ( 392/ 990- 454/1114)  81.9 E40-E36 N19-N26
      9 1731  944 (S18W43 L= 63)   844    201        2 (1644/ 919-1835/ 970) -21.3 W35-W53 S22-S14
     10 1357  890 (S10W15 L= 36)   361     95        1 (1326/ 843-1383/ 934)  38.7 W12-W17 S11-S07
     11  477  864 (N11E37 L=343)   296    100        1 ( 457/ 828- 504/ 916) -80.5 E38-E36 N08-N14
     12  223  778 (N10E64 L=317)   629    100        1 ( 203/ 740- 252/ 830) -63.4 E65-E60 N07-N13
     13 1234  650 (S20W02 L= 22)   352     83        1 (1215/ 612-1258/ 697)  80.5 W00-W03 S22-S16
     14  880  583 (S15E20 L=360)   639    218        0 ( 843/ 490- 929/ 673)  45.0 E25-E15 S19-S09
     15  746  497 (S17E32 L=349)   250    106        0 ( 697/ 484- 799/ 514) -33.7 E34-E28 S18-S14
     16  647  422 (S19E42 L=338)   448    117        0 ( 630/ 375- 678/ 469)  56.3 E44-E38 S22-S16
     17  945  383 (S29E25 L=356)   619     86        1 ( 909/ 350- 989/ 415)  14.0 E28-E20 S30-S27
     18  586  356 (S22E51 L=329)   875     77        2 ( 548/ 325- 618/ 393) -53.1 E55-E47 S24-S19
  

    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