Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1206 1609 (N32E06 L=139)  1409    342        1 (1136/1502-1297/1671)  90.0 E09-W01 N23-N37
      2  553 1603 (N15E47 L= 98)   794    174        2 ( 533/1511- 582/1669) -45.0 E52-E42 N09-N20
      3 1530 1566 (N41W19 L=164)   906    252        0 (1432/1476-1610/1636) -16.4 W09-W26 N37-N43
      4  683 1590 (N17E37 L=108)   417    144        0 ( 615/1576- 760/1598)  24.0 E41-E32 N15-N19
      5  726 1497 (N12E30 L=115)   305    129        0 ( 713/1434- 742/1556) -54.5 E32-E27 N08-N16
      6 1193  704 (S20W19 L=164)   486    197        1 (1144/ 610-1227/ 773) -45.0 W13-W23 S25-S17
      7  491  714 (S39E28 L=117)   234     51        1 ( 462/ 702- 521/ 730)   9.5 E31-E25 S39-S37
      8  743  618 (S38E06 L=139)   279     96        1 ( 720/ 574- 778/ 664) -33.7 E08-E02 S40-S35
      9 1514  577 (S17W44 L=189)   297    127        0 (1461/ 538-1545/ 597)   0.0 W40-W48 S19-S16
     10 1553  464 (S22W54 L=199)   389     91        0 (1515/ 445-1594/ 494)  53.1 W51-W56 S23-S18
     11 1544  315 (S29W69 L=214)  1259    162        2 (1473/ 296-1628/ 334)  19.4 W60-W78 S32-S25
     12  826  229 (S65W22 L=167)  1284    131        4 ( 784/ 187- 875/ 287) -14.7 W13-W32 S68-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