Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  474 1606 (N39E49 L=237)   573    125        1 ( 448/1571- 502/1633) -26.6 E52-E44 N36-N41
      2  903 1500 (N33E10 L=275)   255     47        2 ( 878/1479- 936/1518)  26.6 E11-E07 N31-N34
      3  210 1316 (N18E65 L=220)   587    104        3 ( 191/1292- 231/1337)  56.3 E68-E62 N16-N19
      4  742 1211 (N14E19 L=267)   755    245        2 ( 657/1141- 809/1303)  45.0 E23-E14 N09-N19
      5  178  962 (S04E65 L=220)   721    123        3 ( 155/ 929- 207/1004) -33.7 E68-E61 S06-S01
      6  650  874 (S07E23 L=262)   422    137        1 ( 599/ 836- 698/ 901)   9.5 E26-E20 S09-S05
      7  208  819 (S13E62 L=223)  1156    264        3 ( 171/ 774- 255/ 880)  37.9 E67-E55 S16-S08
      8 1855  818 (S10W65 L=350)   746    123        1 (1820/ 783-1886/ 847)  18.4 W60-W69 S11-S07
      9  686  705 (S18E22 L=264)   243     40        0 ( 659/ 684- 709/ 727) -33.7 E23-E19 S19-S16
     10  938  609 (S24E05 L=280)  1164    302        1 ( 826/ 554- 999/ 721)  47.7 E12-W00 S27-S15
     11  432  609 (S26E42 L=243)   319    101        2 ( 417/ 552- 448/ 643)  80.5 E43-E40 S29-S23
     12  671  572 (S27E24 L=261)   302    112        0 ( 653/ 528- 687/ 617)  74.1 E25-E22 S30-S23
     13  420  393 (S43E53 L=232)  1245    262        1 ( 372/ 370- 481/ 411)  14.9 E60-E45 S45-S40
  

    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