Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1530 1759 (N43W60 L=123)   873    120        0 (1500/1718-1561/1796)  71.6 W59-W61 N39-N45
      2  539 1356 (N33E30 L= 32)   365    138        2 ( 489/1338- 591/1384) -20.6 E34-E26 N30-N34
      3  264 1277 (N29E54 L=  8)  2195    355        5 ( 174/1219- 355/1343)  20.6 E64-E45 N25-N33
      4 1589 1295 (N13W42 L=105)   457    186        1 (1512/1266-1632/1343)  10.3 W35-W46 N12-N15
      5  703 1086 (N14E19 L= 43)   497    136        1 ( 645/1040- 762/1116) -40.6 E23-E15 N10-N16
      6  221  690 (S07E67 L=356)  2610    358        5 ( 179/ 569- 256/ 820)  79.4 E72-E62 S15-N01
      7  345  674 (S09E53 L= 10)  1441    291        0 ( 310/ 568- 420/ 746)  59.7 E58-E45 S15-S03
      8 1787  637 (S32W53 L=116)   531    125        0 (1755/ 615-1828/ 664)   8.1 W50-W57 S33-S30
      9 1222  604 (S23W07 L= 70)   394     86        1 (1173/ 575-1273/ 634)   8.1 W03-W10 S24-S21
  

    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