Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1351 1753 (N43W48 L= 82)   482    112        1 (1328/1716-1378/1776)  12.5 W43-W52 N41-N44
      2 1238 1402 (N25W22 L= 56)   426    123        2 (1192/1353-1282/1439)  18.4 W17-W26 N22-N26
      3  863 1380 (N31E03 L= 31)   316    106        1 ( 808/1366- 920/1395) -18.4 E06-W00 N29-N32
      4  357 1349 (N38E44 L=350)  1025    207        1 ( 297/1302- 427/1384) -17.1 E50-E37 N34-N41
      5  408 1095 (N22E40 L=354)   721    279        1 ( 294/1054- 497/1120) -17.4 E49-E33 N17-N24
      6  688  987 (N11E21 L= 13)   332     65        1 ( 651/ 970- 725/1009) -21.8 E23-E18 N09-N12
      7  223  926 (N13E59 L=335)   596    138        2 ( 186/ 889- 258/ 980) -51.3 E62-E55 N10-N16
      8 1422  385 (S44W15 L= 49)   571    114        1 (1368/ 370-1471/ 411) -18.4 W10-W19 S47-S41
  

    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