Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1848 (N57E07 L= 33)  1272    125        1 ( 865/1808- 969/1891)  39.3 E12-E01 N52-N62
      2 1795 1435 (N22W62 L=102)  4266    246        6 (1698/1347-1889/1497) -15.3 W52-W76 N16-N25
      3  628 1356 (N19E25 L= 15)   398    155        0 ( 581/1291- 654/1395) -66.8 E28-E22 N14-N21
      4  187 1040 (N02E61 L=339)   783    113        0 ( 133/ 997- 230/1080)  18.4 E68-E55 N00-N04
      5  217  657 (S21E67 L=333)  2990    241        4 ( 170/ 557- 270/ 744) -76.0 E74-E61 S27-S15
      6 1820  689 (S24W63 L=103)   666     94        0 (1793/ 648-1843/ 741)  80.5 W60-W64 S26-S20
      7  923  630 (S27E08 L= 32)   227     83        0 ( 882/ 617- 972/ 643)  -9.5 E10-E04 S27-S25
      8 1377  518 (S35W25 L= 65)  1001    279        2 (1317/ 403-1448/ 630)  72.6 W22-W28 S43-S27
      9  621  576 (S29E30 L=  9)   513    142        0 ( 552/ 547- 689/ 607)   9.5 E36-E24 S30-S27
     10 1005  341 (S48E05 L= 35)   329     73        0 ( 969/ 323-1040/ 361) -26.6 E07-E01 S49-S46
     11  771  287 (S52E29 L= 11)   689    116        0 ( 721/ 259- 815/ 324)  21.0 E35-E22 S53-S48
  

    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