Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1625 (N31W49 L= 31)   876    192        1 (1401/1600-1556/1645) -24.4 W42-W54 N27-N33
      2  823 1502 (N40E02 L=340)   621    110        0 ( 775/1449- 876/1544)  20.6 E06-W02 N37-N42
      3 1607 1346 (N12W45 L= 27)   438     86        0 (1585/1302-1628/1381) -90.0 W43-W47 N08-N13
      4  807 1335 (N30E07 L=335)   319     94        2 ( 767/1311- 863/1356) -26.6 E10-E03 N27-N31
      5 1225 1305 (N19W18 L=  0)  1007    351        3 (1071/1278-1398/1322) -25.3 W08-W28 N13-N22
      6 1033 1115 (N12W02 L=344)   309    131        0 ( 981/1057-1074/1159) -66.8 W00-W03 N08-N15
      7 1187  631 (S20W02 L=343)   339    145        0 (1147/ 550-1220/ 679)  49.4 E02-W04 S24-S16
      8 1006  362 (S34E18 L=324)   410     60        2 ( 974/ 339-1034/ 387) -45.0 E20-E15 S36-S32
      9 1390  290 (S53W10 L=352)   397     57        2 (1356/ 275-1420/ 303) -15.9 W06-W13 S55-S51
  

    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