Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  376 1578 (N11E60 L=207)   581     60        1 ( 350/1558- 407/1600)  -8.1 E63-E56 N09-N12
      2 1708 1409 (N35W37 L=303)  1472    261        1 (1594/1341-1805/1478)  -2.6 W25-W47 N33-N36
      3  337 1332 (S04E52 L=215)   481    124        0 ( 303/1289- 383/1379)  60.3 E53-E49 S07-N00
      4  339 1236 (S10E49 L=218)   370     60        0 ( 322/1206- 360/1270) -45.0 E50-E46 S11-S07
      5  965 1177 (S00E07 L=260)   188     51        0 ( 952/1150- 978/1215) -33.7 E08-E05 S01-N01
      6 1874 1052 (N19W57 L=324)   478     29        1 (1850/1029-1898/1077)  53.1 W54-W60 N16-N20
      7 1223  781 (S15W18 L=284)   262    113        0 (1197/ 724-1261/ 821)  90.0 W16-W18 S18-S12
      8 1204  664 (S22W20 L=287)   470    129        0 (1176/ 621-1256/ 721) -20.6 W16-W24 S23-S19
      9 1186  594 (S27W22 L=288)   215     58        0 (1164/ 566-1205/ 624) -26.6 W19-W24 S28-S25
     10  734  487 (S46E04 L=262)  1901    428        3 ( 566/ 417- 934/ 563)   3.5 E20-W13 S48-S43
  

    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