Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1017 1868 (N48E38 L=220)   744     89        0 ( 985/1834-1058/1888)  53.1 E41-E34 N44-N51
      2  925 1811 (N39E38 L=220)   767     98        0 ( 905/1756- 953/1857) -59.0 E42-E33 N34-N44
      3 1265 1552 (N32E02 L=256)   681     84        2 (1244/1495-1283/1608) -69.4 E03-W00 N28-N36
      4 1500 1472 (N35W17 L=274)  1526    382        0 (1397/1309-1596/1597) -31.0 W06-W27 N26-N41
      5 1662 1061 (N15W37 L=295)   872    126        3 (1614/1013-1709/1112) -20.6 W33-W41 N12-N17
      6 1589  852 (N02W38 L=296)   621    110        2 (1540/ 804-1636/ 898)  63.4 W35-W40 S01-N05
      7  919  615 (S32W07 L=264)   329     73        1 ( 891/ 584- 937/ 660) -80.5 W05-W08 S34-S28
      8  719  619 (S38E07 L=251)   289     60        1 ( 694/ 596- 749/ 644)  53.1 E08-E05 S39-S35
      9  888  483 (S41W10 L=268)   264     60        0 ( 859/ 458- 912/ 512)  90.0 W08-W11 S43-S38
     10  772  446 (S48W03 L=261)   302     67        0 ( 734/ 431- 809/ 464)  18.4 W00-W06 S48-S46
     11 1158  339 (S40W37 L=294)   428     74        0 (1121/ 317-1196/ 365)   8.1 W32-W40 S41-S39
     12 1413  331 (S31W55 L=313)   864     97        0 (1363/ 304-1457/ 357)   5.7 W49-W59 S32-S29
  

    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