Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  995 1621 (N34E17 L=180)  2696    614        9 ( 747/1500-1225/1744)  -5.7 E39-W01 N30-N38
      2  397 1595 (N20E60 L=136)   896    328        2 ( 380/1501- 417/1679) -45.0 E67-E54 N14-N25
      3 1299 1496 (N33W08 L=204)   275     58        0 (1270/1478-1333/1512)  21.8 W05-W10 N31-N34
      4  621 1445 (N15E35 L=162)   353     79        0 ( 596/1409- 650/1487)  90.0 E35-E33 N12-N18
      5  385 1270 (N00E47 L=150)   460     87        1 ( 346/1244- 422/1305)  36.9 E49-E44 S01-N02
      6 1598 1074 (N13W35 L=231)   274     54        2 (1576/1041-1616/1110) -63.4 W33-W36 N11-N15
      7 1867 1001 (N15W61 L=258)   740    172        1 (1823/ 969-1899/1035) -10.3 W55-W66 N12-N16
      8  766  959 (S11E14 L=182)   230     72        1 ( 740/ 931- 797/ 988) -21.8 E16-E11 S11-S09
      9  931  909 (S10E03 L=193)   333     82        2 ( 889/ 881- 979/ 934)  31.0 E05-W00 S12-S07
     10 1674  681 (S08W51 L=248)   594    134        2 (1661/ 626-1691/ 746) -60.3 W48-W54 S11-S04
     11 1340  623 (S19W29 L=225)   235     74        2 (1309/ 610-1381/ 638)  36.9 W26-W30 S19-S16
     12 1010  569 (S29W10 L=206)   269     47        0 ( 988/ 547-1035/ 596) -71.6 W07-W11 S30-S27
  

    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