Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1522 1585 (N21W51 L=208)  1056    184        1 (1454/1532-1607/1621) -12.1 W44-W58 N18-N22
      2 1127 1566 (N33W24 L=181)   564     99        0 (1074/1540-1173/1590)  -6.3 W19-W28 N31-N34
      3 1753 1345 (N02W57 L=213)  1000    172        0 (1719/1260-1784/1419)  58.0 W54-W60 S03-N05
      4  674 1393 (N36E12 L=145)   212     48        0 ( 640/1376- 702/1403) -11.3 E14-E09 N35-N37
      5  151  930 (N20E68 L= 88)  1008    178        2 ( 127/ 844- 176/ 997) -90.0 E73-E63 N15-N24
      6 1277  743 (S16W07 L=163)   381    146        1 (1226/ 698-1323/ 805)  20.6 W02-W10 S18-S14
      7  831  693 (S07E20 L=137)  1067    261        2 ( 756/ 581- 910/ 798) -86.2 E21-E17 S15-N00
      8 1329  556 (S30W05 L=161)   385     99        0 (1280/ 541-1378/ 576) -15.9 W01-W08 S31-S27
      9  974  362 (S32E24 L=132)  1037    192        2 ( 914/ 279-1014/ 425)  72.9 E28-E21 S38-S25
  

    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