Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  554 1746 (N29E57 L=142)  1306    172        3 ( 515/1710- 604/1774)  38.7 E61-E51 N25-N31
      2  366 1572 (N13E63 L=136)  1337    159        4 ( 330/1517- 397/1630) -63.4 E67-E59 N09-N17
      3  826 1264 (N05E17 L=182)   311    150        0 ( 770/1203- 855/1330)  90.0 E18-E16 N00-N09
      4 1861  862 (N10W63 L=262)   697    143        2 (1840/ 834-1884/ 901) -14.0 W59-W67 N08-N11
      5 1165  832 (S11W13 L=212)   939    253        0 (1074/ 765-1266/ 880)  45.0 W09-W19 S16-S06
      6  738  824 (S22E12 L=187)   342    154        0 ( 698/ 781- 806/ 848)  39.8 E14-E08 S24-S19
      7 1479  796 (S05W33 L=232)   468     96        1 (1449/ 752-1499/ 846) -71.6 W30-W34 S07-S01
      8  289  654 (S45E48 L=151)   454     99        1 ( 255/ 630- 326/ 680)   0.0 E51-E43 S45-S43
      9  608  619 (S38E17 L=182)  1079    319        0 ( 484/ 592- 762/ 641)  18.4 E27-E06 S41-S33
     10  296  561 (S52E50 L=149)   525    242        0 ( 249/ 543- 336/ 582)  21.0 E57-E44 S54-S49
     11 1388  453 (S26W41 L=240)  2387    704        2 (1061/ 382-1615/ 497)  14.0 W18-W62 S34-S22
  

    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