Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1006 1632 (N37W02 L=324)   352     67        1 ( 985/1598-1021/1669)  68.2 W00-W02 N34-N39
      2 1635 1403 (N19W43 L=  6)   322     82        2 (1610/1367-1663/1438) -78.7 W42-W45 N16-N21
      3 1736 1376 (N17W52 L= 14)   380     78        1 (1706/1354-1768/1397) -31.0 W49-W54 N15-N18
      4 1464 1189 (N06W28 L=350)  1384    305        2 (1337/1048-1541/1279) -49.8 W20-W32 S02-N11
      5 1795 1150 (N03W54 L= 16)   346     71        1 (1773/1112-1809/1187)  78.7 W51-W55 N00-N05
      6 1009 1160 (N06E00 L=322)   231     84        2 ( 974/1120-1042/1184) -71.6 E02-W01 N03-N07
      7  418  790 (S14E41 L=281)  2154    560        7 ( 331/ 622- 547/ 956)  42.1 E51-E30 S23-S04
      8 1129  625 (S27W05 L=328)   207     65        0 (1104/ 597-1152/ 664)  59.0 W03-W06 S29-S24
      9 1461  403 (S44W34 L=356)   372     82        1 (1430/ 388-1500/ 420)  -8.1 W30-W37 S45-S42
  

    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