Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1026 1624 (N36W05 L= 32)   261     70        1 (1009/1591-1049/1652)  78.7 W03-W06 N33-N38
      2 1748 1392 (N16W54 L= 82)  1693    357        7 (1701/1262-1802/1472)  77.0 W48-W61 N08-N21
      3 1223 1219 (N08W13 L= 41)   234     73        0 (1180/1206-1257/1236) -21.8 W10-W15 N07-N09
      4  253 1086 (N07E53 L=335)   557    125        2 ( 220/1036- 289/1149)  54.5 E56-E49 N04-N11
      5  425 1044 (N03E38 L=350)  1161    266        3 ( 367/ 913- 463/1143) -72.9 E42-E35 S04-N09
      6 1882 1044 (S05W63 L= 91)   423     82        1 (1866/1001-1902/1074)  45.0 W61-W65 S07-S03
      7  197  989 (N02E59 L=329)  2602    443       11 ( 124/ 919- 289/1049) -18.4 E69-E49 S02-N06
      8  275  770 (S11E54 L=333)   899    208        1 ( 228/ 712- 323/ 830)  28.6 E60-E48 S14-S08
      9  362  605 (S22E51 L=337)   537     85        0 ( 338/ 556- 394/ 648) -80.5 E52-E49 S25-S19
     10  973  401 (S42E09 L= 18)   814    204        1 ( 888/ 361-1073/ 459)  20.6 E16-W00 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