Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  663 1148 (N20E21 L=196)   530    192        1 ( 582/1117- 747/1176) -15.3 E26-E14 N17-N22
      2 1575  954 (S08W33 L=250)   533     91        3 (1539/ 910-1610/1000) -66.8 W31-W35 S11-S04
      3  778  694 (S09E21 L=196)   521    105        1 ( 739/ 658- 814/ 739)  59.0 E23-E18 S11-S06
      4  305  601 (S10E61 L=156)   874    121        1 ( 276/ 554- 344/ 660) -76.0 E63-E58 S13-S05
      5  350  627 (S08E55 L=162)   503     71        2 ( 328/ 594- 378/ 649)  23.2 E58-E51 S10-S06
      6 1672  565 (S36W41 L=259)  1469    186        3 (1585/ 525-1742/ 613)   0.0 W32-W49 S38-S34
      7  543  553 (S15E41 L=176)   271     92        1 ( 517/ 512- 566/ 581)  76.0 E42-E39 S17-S13
      8 1405  515 (S33W17 L=234)   739    148        0 (1342/ 484-1459/ 556)  11.3 W11-W21 S34-S31
      9  495  503 (S18E47 L=170)   646    128        1 ( 448/ 472- 545/ 548) -51.3 E49-E42 S20-S14
     10 1024  448 (S30E11 L=206)  1042    243        1 ( 948/ 376-1114/ 517)  25.0 E18-E03 S34-S25
     11  704  484 (S22E32 L=186)   437     89        1 ( 674/ 456- 754/ 511)   0.0 E33-E27 S23-S19
     12  557  330 (S32E53 L=164)  2084    232        1 ( 526/ 239- 605/ 410)  26.6 E73-E43 S40-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