Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  995 1544 (N31W04 L=109)   274     70        2 ( 968/1517-1022/1584) -90.0 W01-W05 N29-N34
      2  265 1339 (N26E52 L= 53)  2162    403        2 ( 202/1188- 307/1447) -63.4 E60-E47 N15-N33
      3 1338 1251 (N09W21 L=126)   362     88        1 (1293/1207-1374/1283) -51.3 W18-W23 N06-N11
      4  212 1203 (N18E56 L= 49)   369     83        1 ( 197/1165- 229/1245)  59.0 E57-E54 N15-N20
      5 1785  749 (S24W52 L=157)   360    110        0 (1746/ 737-1821/ 766)   0.0 W48-W55 S25-S23
      6  489  724 (S14E38 L= 67)   295     75        2 ( 461/ 706- 533/ 747) -21.8 E39-E34 S15-S12
      7 1636  532 (S38W42 L=147)   541    100        1 (1587/ 491-1685/ 570)  35.5 W38-W46 S40-S35
      8  407  481 (S28E52 L= 53)  1729    286        3 ( 351/ 418- 459/ 563)  34.7 E60-E47 S31-S22
  

    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