Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1178 1855 (N65E07 L= 81)  4914    418        7 ( 992/1770-1350/1928)  -4.2 E43-W17 N58-N71
      2 1120 1685 (N45E06 L= 83)  1644    307        4 ( 972/1628-1220/1724)  -2.6 E19-W03 N42-N47
      3  555 1668 (N33E52 L= 37)  1286    264        2 ( 482/1583- 663/1708)  33.7 E57-E42 N26-N37
      4  347 1402 (N13E55 L= 33)  1684    357        7 ( 298/1277- 399/1542) -46.8 E64-E48 N05-N21
      5 1534 1188 (N18W31 L=119)  1113    299        0 (1444/1074-1612/1297) -40.6 W23-W37 N11-N23
      6 1491 1150 (N15W28 L=116)   228     87        0 (1453/1123-1531/1170)  38.7 W25-W30 N12-N16
      7 1861 1040 (N13W63 L=151)  1031    198        1 (1837/ 956-1891/1127)  84.3 W60-W68 N08-N18
      8  380  994 (S11E43 L= 46)   418    134        1 ( 328/ 945- 432/1035)  36.9 E46-E38 S15-S08
      9  160  859 (S23E69 L= 19)   847     72        1 ( 143/ 823- 181/ 895) -90.0 E74-E65 S25-S20
     10  274  782 (S26E53 L= 35)   444     84        0 ( 254/ 739- 298/ 825)  45.0 E56-E50 S29-S23
     11  759  751 (S21E13 L= 75)   389     94        1 ( 716/ 727- 804/ 773)  -9.5 E15-E09 S21-S19
     12 1137  649 (S21W13 L=101)   733    332        1 (1018/ 548-1256/ 749)  55.5 W07-W18 S29-S13
     13  585  708 (S26E25 L= 63)   833    175        2 ( 501/ 685- 670/ 737)  17.1 E31-E18 S28-S23
     14 1230  602 (S23W20 L=109)   297    150        0 (1190/ 546-1279/ 660) -33.7 W16-W25 S25-S19
     15  684  212 (S69E25 L= 64)   884     79        3 ( 657/ 176- 714/ 254)  76.0 E31-E19 S75-S63
     16  946  180 (S64W18 L=106)   591     84        2 ( 916/ 154- 986/ 212) -12.1 W12-W26 S67-S60
  

    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