Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1279 1627 (N44W09 L=318)   658    129        1 (1238/1572-1308/1674) -53.1 W04-W11 N39-N47
      2  365 1602 (N26E68 L=241)   637     72        1 ( 346/1571- 390/1640) -59.0 E71-E64 N23-N28
      3 1781 1345 (N32W56 L=  5)  1020    227        0 (1717/1252-1843/1424) -39.8 W50-W63 N26-N36
      4 1770 1256 (N26W53 L=  2)   414     65        1 (1746/1226-1802/1278) -18.4 W50-W56 N24-N26
      5 1272 1208 (N15W13 L=321)   335     90        0 (1233/1175-1313/1249) -26.6 W09-W15 N13-N16
      6 1869 1210 (N24W66 L= 15)   848     58        1 (1842/1188-1892/1247)  26.6 W61-W70 N22-N26
      7 1629  662 (S13W48 L=357)   807    235        0 (1557/ 598-1685/ 765)  67.4 W44-W50 S17-S05
      8 1414  675 (S15W31 L=340)   879    229        2 (1348/ 601-1457/ 733) -77.5 W25-W34 S20-S11
      9 1500  640 (S16W38 L=347)   251     52        0 (1485/ 611-1514/ 664) -71.6 W36-W39 S17-S14
     10  779  592 (S31E10 L=299)  1477    424        3 ( 611/ 540- 981/ 648)  25.6 E22-W03 S37-S23
     11 1590  379 (S32W63 L= 12)  2059    195        0 (1509/ 354-1680/ 426)  19.4 W53-W74 S34-S27
     12 1122  333 (S44W23 L=332)  1299    203        3 (1038/ 307-1227/ 360)  23.2 W15-W30 S47-S40
  

    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