Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1333 1679 (N35W46 L=149)  3298    562        3 (1192/1508-1542/1753)  -7.1 W28-W69 N27-N39
      2  571 1135 (N24E25 L= 79)   390     85        3 ( 540/1099- 601/1175) -71.6 E27-E23 N21-N27
      3 1732  898 (S21W42 L=146)  1146    291        1 (1626/ 850-1800/ 992)  -3.4 W33-W50 S23-S15
      4 1509  944 (S11W26 L=130)   227     51        1 (1482/ 919-1533/ 970) -63.4 W24-W27 S12-S08
      5  227  875 (N16E60 L= 44)   361     58        1 ( 210/ 844- 242/ 905)  53.1 E62-E57 N13-N17
      6 1854  849 (S28W56 L=159)  1479    155        2 (1788/ 812-1919/ 900)   0.0 W46-W67 S30-S26
      7  170  799 (N11E71 L= 33)  2631    152        4 ( 122/ 733- 213/ 850) -28.3 E79-E65 N07-N14
      8 1407  752 (S19W15 L=119)   498    168        0 (1391/ 678-1421/ 836)  60.9 W12-W17 S22-S13
      9 1760  697 (S35W44 L=147)  1101    130        4 (1729/ 640-1793/ 762)  45.0 W40-W47 S38-S30
     10  372  566 (S05E57 L= 47)   486     74        1 ( 331/ 547- 405/ 585)   0.0 E60-E52 S06-S03
     11 1260  476 (S32E01 L=103)   270     72        0 (1235/ 440-1276/ 515)  53.1 E03-W01 S34-S29
     12  745  349 (S27E39 L= 64)   373     47        1 ( 715/ 326- 774/ 373)  18.4 E42-E36 S27-S25
  

    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