Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1820 (N62W27 L=199)  1784    307        3 ( 758/1747- 988/1885)   0.0 W09-W48 N58-N65
      2  596 1655 (N60E15 L=157)   658     89        0 ( 575/1611- 613/1706)  69.4 E17-E11 N55-N64
      3  665 1080 (N19E20 L=152)   337    153        0 ( 626/1031- 730/1109)   8.1 E22-E15 N16-N19
      4 1605  877 (S18W31 L=203)   351     65        2 (1576/ 852-1630/ 903)  11.3 W28-W33 S19-S16
      5 1678  855 (S22W36 L=208)   346     91        0 (1640/ 833-1700/ 892)   0.0 W33-W38 S22-S19
      6 1460  751 (S21W18 L=191)   338     98        0 (1442/ 715-1477/ 799)  53.1 W16-W20 S23-S18
      7  811  582 (S13E24 L=148)   392     82        0 ( 777/ 567- 860/ 603) -36.9 E26-E21 S14-S10
      8 1036  544 (S21E13 L=159)   261    114        0 (1003/ 500-1071/ 599)  26.6 E16-E08 S23-S18
      9  694  528 (S13E34 L=138)   686    129        3 ( 626/ 491- 757/ 567)   0.0 E39-E28 S15-S11
     10 1269  421 (S37E03 L=169)   404     94        1 (1228/ 402-1321/ 441)  -9.5 E05-W01 S37-S34
  

    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