Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  430 1600 (N38E53 L=238)   590    131        1 ( 403/1574- 468/1629) -12.5 E56-E47 N36-N40
      2  198 1292 (N16E66 L=224)   887    200        0 ( 167/1264- 233/1332) -23.2 E72-E62 N14-N18
      3  664 1215 (N14E24 L=267)   679    229        1 ( 580/1150- 729/1300)  45.0 E29-E19 N09-N19
      4  532 1184 (N11E33 L=258)   357     79        2 ( 508/1159- 568/1212) -26.6 E34-E30 N09-N13
      5  158  912 (S08E68 L=222)  1439    128        1 ( 135/ 849- 175/ 967) -66.8 E72-E65 S11-S04
      6 1777  769 (S12W56 L=347)  3251    493        3 (1654/ 559-1886/ 962)  54.2 W48-W68 S26-N00
      7  559  889 (S07E30 L=261)   317    106        0 ( 519/ 861- 608/ 906) -18.4 E32-E26 S08-S05
      8  279  822 (S12E54 L=237)   459     83        1 ( 257/ 779- 302/ 861) -68.2 E55-E51 S15-S09
      9  873  646 (S21E09 L=282)  1278    365        2 ( 754/ 555- 961/ 786)  49.1 E17-E03 S27-S11
     10  611  709 (S18E27 L=264)   229     36        0 ( 591/ 689- 630/ 733) -56.3 E28-E25 S19-S16
     11  599  585 (S26E29 L=262)   282    117        0 ( 582/ 535- 615/ 634)  63.4 E31-E27 S30-S22
     12  398  415 (S41E54 L=237)  1265    297        2 ( 352/ 386- 446/ 450)  23.6 E62-E46 S44-S37
  

    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