Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1026 1859 (N51E23 L=293)   655     88        1 ( 989/1835-1072/1885)  -5.7 E29-E19 N48-N52
      2  755 1698 (N31E33 L=283)  1233    353        1 ( 705/1548- 818/1819) -67.8 E39-E29 N19-N41
      3  969 1690 (N34E18 L=298)   373     98        0 ( 949/1638- 983/1738) -49.4 E21-E15 N30-N37
      4 1512 1508 (N33W24 L=340)   272     89        0 (1469/1496-1557/1519)  23.2 W20-W27 N31-N34
      5  623  943 (S18E23 L=293)   633    180        1 ( 569/ 869- 684/ 992)  53.1 E26-E19 S22-S14
      6  423  578 (S44E40 L=276)  1143    174        4 ( 376/ 529- 472/ 637) -60.3 E45-E34 S47-S40
      7 1738  472 (S22W72 L= 28)  1269    105        2 (1712/ 426-1766/ 514) -26.6 W65-W78 S24-S20
      8 1096  216 (S57W36 L=353)   969    142        1 (1027/ 180-1154/ 248)  13.0 W26-W45 S59-S54
  

    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