Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1573 1654 (N32W61 L=325)  2219    193        4 (1498/1616-1653/1691) -23.2 W51-W70 N27-N34
      2  696 1497 (N42E15 L=248)  1146    261        1 ( 569/1477- 812/1514)  -8.1 E26-E05 N39-N44
      3  276 1460 (N46E63 L=200)  2132    151        2 ( 203/1428- 354/1484)   1.7 E84-E50 N43-N46
      4 1242 1019 (N03W13 L=276)   205    116        0 (1210/ 985-1288/1063) -59.0 W11-W14 N00-N05
      5  207  877 (N09E62 L=201)   668    132        1 ( 189/ 815- 226/ 945) -82.9 E63-E60 N05-N13
      6 1886  764 (S29W66 L=329)   805     76        1 (1863/ 733-1904/ 803)  36.9 W61-W70 S31-S26
      7 1811  656 (S34W56 L=319)  1754    131        4 (1743/ 614-1864/ 692) -11.3 W46-W66 S37-S30
      8  310  562 (S11E62 L=201)  3183    254        6 ( 256/ 431- 345/ 660)  38.7 E76-E54 S21-S05
      9  772  513 (S20E26 L=237)   265     61        0 ( 754/ 486- 789/ 534) -76.0 E27-E24 S22-S18
     10 1174  309 (S44E05 L=258)  3140    525        4 ( 980/ 227-1397/ 416)   6.2 E23-W14 S47-S39
     11  563  330 (S30E52 L=211)   953    167        1 ( 512/ 267- 593/ 385)  63.4 E58-E47 S35-S25
     12  892  200 (S47E35 L=229)   546     63        0 ( 862/ 181- 925/ 223) -53.1 E37-E32 S49-S44
  

    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