Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  457 1688 (N19E63 L= 44)   566     84        2 ( 429/1645- 475/1720) -24.4 E69-E57 N16-N21
      2 1048 1421 (N16E09 L= 97)   466    171        0 (1006/1345-1080/1496) -41.2 E13-E05 N12-N19
      3 1179 1372 (N16E01 L=106)   364    105        0 (1135/1356-1225/1397)  26.6 E03-W01 N15-N18
      4 1283 1279 (N14W08 L=114)   229    107        0 (1266/1230-1309/1329) -51.3 W05-W10 N11-N16
      5 1713  792 (N00W50 L=156)   663     90        2 (1687/ 745-1737/ 835) -26.6 W46-W53 S01-N02
      6 1331  734 (S14W26 L=132)   205     60        1 (1310/ 714-1363/ 759) -11.3 W23-W28 S15-S13
      7  972  532 (S37W12 L=119)  1231    315        4 ( 854/ 495-1108/ 573)  33.7 W04-W19 S41-S30
      8 1464  462 (S25W48 L=154)   317     64        2 (1439/ 434-1479/ 493) -76.0 W45-W49 S26-S22
      9 1153  412 (S38W30 L=137)   279     80        0 (1140/ 377-1169/ 456) -56.3 W28-W32 S40-S34
     10 1085  316 (S46W33 L=140)   350     79        0 (1048/ 305-1138/ 333)  26.6 W29-W36 S48-S44
  

    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