Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  768 1648 (N30E35 L= 98)   326     62        0 ( 742/1619- 787/1669) -21.8 E37-E32 N28-N31
      2  494 1535 (N16E49 L= 83)   575    148        4 ( 441/1504- 530/1562)  39.8 E52-E45 N13-N18
      3  766 1456 (N17E26 L=106)   621    164        3 ( 700/1419- 842/1479)  32.0 E30-E21 N14-N20
      4 1362 1390 (N27W13 L=145)   444     82        1 (1332/1344-1381/1428) -78.7 W10-W14 N24-N29
      5  259 1126 (S12E55 L= 77)   365    130        1 ( 234/1062- 278/1182) -45.0 E59-E52 S15-S09
      6  589  983 (S14E26 L=107)   208     50        0 ( 576/ 954- 604/1011) -56.3 E26-E24 S15-S12
      7  486  843 (S25E31 L=101)  1025    249        2 ( 394/ 773- 553/ 939)  51.8 E38-E26 S31-S17
      8 1344  537 (S23W33 L=166)   588    112        1 (1308/ 502-1378/ 591)  90.0 W31-W36 S26-S19
      9 1574  536 (S18W50 L=183)   303     93        0 (1548/ 502-1595/ 581) -80.5 W49-W52 S20-S14
     10  680  456 (S47E11 L=122)  1428    249        4 ( 571/ 432- 811/ 488)  15.9 E21-W00 S50-S41
     11 1082  195 (S56W41 L=174)  1129    130        0 (1034/ 154-1122/ 248) -50.7 W36-W48 S60-S49
  

    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