Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  841 1763 (N44E21 L= 62)   449     58        0 ( 809/1744- 879/1781)  -9.5 E23-E17 N42-N45
      2  928 1712 (N40E12 L= 71)   258     50        0 ( 914/1686- 944/1738)  90.0 E13-E11 N37-N42
      3  985 1704 (N40E08 L= 75)   362     25        0 ( 960/1687-1013/1723) -14.0 E09-E05 N38-N41
      4  559 1629 (N32E38 L= 44)   980    295        1 ( 481/1550- 621/1698) -32.9 E48-E31 N26-N37
      5 1079 1677 (N38E00 L= 83)   385     46        0 (1056/1650-1110/1698)  71.6 E02-W02 N36-N40
      6 1291 1575 (N32W15 L= 98)   540    101        1 (1234/1553-1346/1599)  -7.1 W11-W19 N30-N33
      7 1502 1517 (N30W31 L=113)   539     92        0 (1462/1496-1545/1534)  15.9 W27-W34 N28-N31
      8 1048 1332 (N14E00 L= 82)   253     78        0 (1035/1292-1060/1379) -80.5 E01-W00 N11-N17
      9 1273 1312 (N14W14 L= 96)   250     70        0 (1233/1278-1302/1333) -26.6 W11-W15 N12-N15
     10 1448  507 (S34W35 L=117)   262     41        0 (1424/ 492-1478/ 524)   0.0 W32-W37 S34-S32
  

    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