Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  660 1147 (N20E21 L=197)   531    210        0 ( 576/1116- 748/1175) -14.0 E26-E14 N17-N22
      2 1576  936 (S09W33 L=251)   630    128        3 (1532/ 883-1618/ 993) -60.3 W30-W35 S12-S04
      3  772  694 (S09E21 L=197)   491     99        1 ( 730/ 661- 806/ 737) -78.7 E24-E19 S11-S06
      4  352  650 (S07E54 L=164)   854    133        3 ( 324/ 590- 376/ 711)  41.2 E58-E49 S10-S03
      5  300  602 (S10E62 L=156)   870    135        0 ( 271/ 554- 341/ 665) -74.1 E64-E59 S13-S05
      6 1665  563 (S36W41 L=258)  1424    193        2 (1574/ 527-1735/ 607)   0.0 W31-W48 S38-S34
      7  509  524 (S17E45 L=173)  1030    255        4 ( 446/ 477- 562/ 583)   0.0 E49-E39 S19-S13
      8 1400  515 (S33W16 L=234)   696    129        1 (1339/ 485-1454/ 550)  12.5 W11-W20 S34-S31
      9 1018  447 (S30E12 L=206)  1047    249        3 ( 941/ 374-1112/ 521)  25.0 E18-E03 S34-S25
     10  694  483 (S22E32 L=185)   370     89        0 ( 668/ 460- 734/ 510)   9.5 E34-E28 S23-S19
     11  554  332 (S31E53 L=164)  2151    247        2 ( 522/ 237- 601/ 428)  27.3 E74-E43 S40-S24
  

    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