Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  477 1652 (N19E56 L=170)   483     48        2 ( 451/1629- 509/1676)  26.6 E58-E52 N16-N20
      2  624 1594 (N17E40 L=185)   404     93        0 ( 601/1561- 656/1639) -26.6 E44-E36 N15-N20
      3  337 1557 (N10E64 L=162)   806     96        2 ( 302/1518- 358/1605) -29.1 E67-E58 N08-N13
      4  509 1489 (N09E44 L=181)   595    150        2 ( 482/1424- 548/1558) -37.9 E47-E38 N05-N12
      5  393 1262 (S07E45 L=180)   328     79        1 ( 350/1241- 430/1282)   0.0 E49-E42 S07-S05
      6  331 1209 (S11E49 L=176)   319     56        0 ( 302/1185- 363/1228)  -9.5 E52-E45 S11-S09
      7  913 1154 (S03E09 L=216)   224     53        1 ( 881/1127- 937/1173) -14.0 E11-E07 S03-S01
      8  380 1048 (S19E42 L=183)   273     28        0 ( 360/1027- 401/1073) -71.6 E43-E40 S20-S17
      9  273  988 (S24E51 L=174)   306     35        0 ( 250/ 969- 291/1007) -14.0 E53-E49 S25-S23
     10  482  948 (S23E32 L=193)   381    123        0 ( 447/ 911- 535/ 998) -15.9 E35-E28 S24-S20
     11  225  884 (S32E57 L=168)   417     73        1 ( 205/ 853- 249/ 915) -38.7 E59-E54 S33-S29
     12  209  786 (S38E61 L=164)   528     64        0 ( 173/ 772- 240/ 805)   0.0 E67-E56 S38-S36
  

    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