Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1102 1862 (N55E10 L=342)   622     96        0 (1067/1836-1161/1877)  18.4 E13-E04 N51-N57
      2  435 1664 (N29E58 L=294)  1147    165        3 ( 402/1593- 464/1740) -47.5 E66-E51 N23-N35
      3 1163 1511 (N26W03 L=355)   407    149        1 (1132/1445-1184/1576) -77.5 W00-W03 N21-N30
      4 1776 1241 (N19W50 L= 42)   954    236        2 (1735/1141-1804/1354) -69.0 W46-W52 N12-N25
      5  381  946 (S17E42 L=310)   305     90        0 ( 367/ 899- 402/ 985)  59.0 E43-E40 S19-S14
      6  612  797 (S24E24 L=328)   441    186        0 ( 578/ 717- 633/ 891)  79.7 E26-E22 S29-S18
      7  199  849 (S24E62 L=290)   773    152        1 ( 168/ 801- 249/ 890)  24.4 E67-E55 S26-S21
      8  436  840 (S23E38 L=314)   252     93        0 ( 398/ 817- 454/ 875)  45.0 E40-E36 S24-S20
      9 1715  653 (S17W55 L= 47)   423     50        0 (1692/ 628-1737/ 676) -36.9 W52-W57 S18-S15
     10  473  596 (S38E37 L=315)   359     86        0 ( 444/ 559- 510/ 641)  36.9 E40-E32 S41-S35
     11 1110  217 (S59W30 L= 22)   435     73        1 (1086/ 193-1136/ 246) -24.0 W27-W36 S61-S56
  

    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