Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1621 (N27W66 L=230)   504     76        1 (1600/1603-1673/1637) -33.7 W63-W69 N25-N29
      2 1397 1465 (N25W35 L=198)  1669    317        4 (1264/1371-1482/1570) -70.6 W30-W39 N17-N34
      3 1261 1171 (N11W17 L=180)   711    319        0 (1201/1107-1351/1216)   5.2 W12-W23 N06-N13
      4 1379  933 (S05W19 L=183)   453    133        1 (1325/ 912-1451/ 956) -35.5 W16-W23 S08-S03
      5 1638  736 (S24W36 L=199)  2226    571        3 (1470/ 564-1785/ 956)  30.4 W21-W50 S31-S13
      6 1619  917 (S12W36 L=199)   350    150        0 (1588/ 886-1650/ 953)  31.0 W33-W38 S13-S09
      7 1488  842 (S13W25 L=188)   395     93        0 (1461/ 795-1507/ 892)  59.0 W22-W26 S15-S10
      8  434  862 (N07E41 L=122)   270     54        1 ( 410/ 843- 460/ 881)   0.0 E43-E38 N05-N07
      9 1738  695 (S30W46 L=209)  1351    240        1 (1669/ 605-1813/ 787)  32.7 W39-W54 S34-S25
     10  615  570 (S14E36 L=128)   342     76        1 ( 578/ 547- 649/ 585)   0.0 E38-E32 S15-S12
     11  544  465 (S20E45 L=118)   408     74        2 ( 526/ 429- 567/ 492)  26.6 E48-E42 S22-S17
  

    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