Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1509 1743 (N31W59 L=202)  2194    437        6 (1413/1691-1612/1776) -10.8 W48-W69 N27-N32
      2 1296 1596 (N29W35 L=178)   345     84        1 (1261/1566-1330/1616)   0.0 W31-W37 N26-N30
      3 1045 1054 (N06W02 L=145)   308    105        0 (1011/1000-1099/1102)  26.6 E01-W05 N03-N07
      4  335  969 (N17E45 L= 98)   272     53        0 ( 320/ 937- 351/1000) -76.0 E46-E43 N15-N19
      5 1841  892 (S26W51 L=195)   691    154        1 (1809/ 823-1874/ 967)  41.2 W47-W56 S29-S22
      6 1604  730 (S28W29 L=172)   762    219        0 (1539/ 630-1661/ 845)  39.8 W22-W34 S32-S22
      7  989  769 (S09E08 L=135)   201     54        0 ( 977/ 736-1004/ 803)  63.4 E09-E06 S10-S06
      8 1173  504 (S30E05 L=138)   326     63        0 (1137/ 487-1203/ 527)   9.5 E08-E02 S30-S28
      9  493  362 (S22E56 L= 87)  2468    500        8 ( 448/ 236- 522/ 499)  48.2 E66-E48 S31-S12
     10  575  438 (S19E45 L= 98)   298     30        0 ( 555/ 416- 595/ 457) -71.6 E46-E43 S20-S16
     11  654  349 (S26E45 L= 99)   690    173        0 ( 635/ 289- 672/ 412)  48.4 E48-E40 S30-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