Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1745 (N54W35 L=211)   471     57        0 (1317/1719-1368/1765)  18.4 W30-W39 N51-N55
      2  841 1732 (N53E19 L=157)   335     42        1 ( 817/1711- 872/1748)  11.3 E21-E16 N50-N54
      3  668 1549 (N37E29 L=148)   654    174        0 ( 593/1494- 731/1580) -23.2 E36-E22 N33-N39
      4 1729 1546 (N36W68 L=245)   696     58        0 (1701/1524-1762/1564)   6.7 W61-W78 N34-N36
      5  832 1494 (N34E14 L=162)   472    122        0 ( 781/1468- 890/1518)  20.6 E18-E10 N31-N35
      6 1043 1456 (N31W01 L=178)   331     81        0 (1003/1432-1081/1475) -36.9 E01-W03 N29-N32
      7  344 1275 (N18E50 L=126)   479    139        1 ( 303/1243- 412/1302)   0.0 E54-E43 N15-N19
      8  387  569 (S27E50 L=126)  1582    358        2 ( 332/ 428- 482/ 715)  32.9 E71-E37 S38-S16
      9 1181  669 (S19W10 L=187)   278    110        1 (1123/ 655-1232/ 690)   0.0 W06-W13 S19-S17
     10  643  493 (S32E29 L=148)  2946    518        8 ( 483/ 348- 817/ 689)  34.0 E53-E13 S44-S17
     11 1492  508 (S31W36 L=213)   323     86        0 (1462/ 483-1531/ 539)  45.0 W33-W38 S32-S28
     12 1413  320 (S46W37 L=214)   485     78        1 (1383/ 295-1429/ 358)  90.0 W34-W40 S48-S42
     13 1465  302 (S48W46 L=223)   550     42        0 (1442/ 283-1489/ 322) -20.0 W40-W51 S50-S46
  

    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