Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1374 1638 (N45W06 L=170)   410    110        0 (1314/1618-1421/1653)  18.4 W01-W10 N42-N46
      2  488 1607 (N17E54 L=109)   542     68        2 ( 459/1580- 525/1636)  -7.1 E58-E50 N15-N19
      3  465 1535 (N12E52 L=112)   435     72        2 ( 430/1505- 490/1567)  76.0 E54-E48 N10-N14
      4 1672 1294 (N30W36 L=200)  2275    400        0 (1529/1195-1810/1423) -12.1 W23-W51 N27-N34
      5  193 1074 (S20E62 L=101)   672    121        0 ( 170/1015- 226/1125)  54.5 E65-E58 S23-S16
      6 1131  999 (S02W07 L=170)   959    261        2 (1078/ 902-1177/1118) -45.0 W00-W10 S07-N02
      7  193  820 (S36E63 L=100)   507     98        0 ( 180/ 784- 211/ 867)  50.2 E66-E60 S38-S32
      8 1060  825 (S14W07 L=171)   288     88        0 (1028/ 803-1102/ 853)  53.1 W05-W08 S15-S11
      9 1462  545 (S19W42 L=206)   406     61        0 (1419/ 526-1495/ 559)   9.5 W38-W45 S20-S17
  

    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