Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  679 1700 (N38E41 L=107)   440     73        1 ( 645/1681- 725/1714)  18.4 E44-E37 N36-N39
      2 1070 1477 (N29E04 L=145)  2194    463        4 ( 872/1369-1262/1568) -13.6 E19-W10 N24-N32
      3 1255 1240 (N16W11 L=160)   248     82        0 (1218/1216-1295/1263)  26.6 W09-W13 N14-N17
      4  292 1194 (N00E53 L= 96)   943    130        2 ( 248/1137- 340/1245)  45.0 E56-E49 S03-N03
      5 1839 1126 (N17W59 L=208)   460     43        0 (1812/1109-1865/1147)  15.9 W55-W62 N16-N18
      6 1310  755 (S12W21 L=170)   630    231        0 (1217/ 708-1398/ 815)  35.0 W16-W26 S15-S07
      7 1398  601 (S20W31 L=180)  1674    524        2 (1189/ 534-1551/ 742)  40.4 W18-W38 S27-S09
      8 1175  258 (S49W31 L=179)   415     56        1 (1136/ 240-1203/ 276)   9.5 W27-W33 S50-S47
  

    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