Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1701 (N51E00 L= 54)   516     89        3 ( 953/1662-1012/1736)  38.7 E02-W03 N47-N54
      2 1124 1603 (N42W12 L= 66)   309     52        0 (1096/1583-1155/1623) -36.9 W09-W13 N40-N43
      3  311 1439 (N32E60 L=355)  1046    141        2 ( 296/1381- 341/1508) -41.6 E66-E53 N28-N36
      4 1413 1487 (N32W32 L= 86)   461    107        0 (1359/1475-1472/1501)  -6.3 W27-W36 N30-N32
      5 1811 1395 (N22W69 L=124)  1235    139        4 (1768/1345-1845/1437) -50.2 W64-W76 N18-N24
      6 1015 1385 (N27W01 L= 56)   313    104        0 ( 968/1364-1063/1401)   9.5 E02-W04 N26-N28
      7  359 1136 (N13E46 L=  8)   472    202        0 ( 311/1067- 427/1207) -39.3 E51-E39 N08-N17
      8 1709  791 (S14W48 L=102)  2652    525        3 (1586/ 610-1796/ 989)  56.9 W38-W55 S25-S02
      9 1024  777 (S11E01 L= 53)   194     53        1 ( 995/ 756-1053/ 797)  26.6 E02-W00 S11-S09
     10 1342  686 (S18W19 L= 74)   250     57        0 (1324/ 661-1362/ 713)  71.6 W18-W20 S20-S16
     11 1103  542 (S27W03 L= 58)   505    173        2 (1036/ 520-1200/ 568)  11.3 E01-W09 S28-S25
     12 1190  164 (S65W16 L= 70)   757     86        0 (1158/ 144-1219/ 185) -22.6 W09-W22 S69-S60
  

    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