Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1221 1630 (N36E07 L= 83)  1925    475        4 (1152/1422-1320/1825) -48.4 E19-W06 N25-N52
      2 1013 1543 (N24E16 L= 73)  1507    469        1 ( 857/1390-1167/1651) -12.5 E28-E01 N18-N28
      3  570 1578 (N15E45 L= 44)   307     86        1 ( 534/1543- 596/1611) -31.0 E47-E41 N13-N16
      4 1447 1403 (N28W15 L=104)   347    106        0 (1430/1361-1480/1451) -51.3 W13-W18 N26-N31
      5 1561 1306 (N26W25 L=114)   817    208        1 (1487/1261-1603/1381)  68.2 W19-W27 N22-N32
      6 1837 1318 (N37W53 L=142)  1063    111        1 (1805/1266-1882/1348)  32.7 W48-W62 N32-N41
      7  263 1311 (S07E59 L= 30)   406    104        0 ( 237/1256- 279/1342) -35.5 E63-E55 S10-S05
      8  704 1264 (S00E24 L= 65)   645    199        2 ( 650/1211- 772/1314)  56.3 E26-E20 S04-N03
      9  454 1278 (S05E41 L= 48)   246     72        1 ( 427/1260- 485/1309)  53.1 E42-E39 S06-S02
     10 1482 1232 (N19W21 L=110)   442    104        0 (1430/1195-1515/1251)  -9.5 W17-W23 N17-N20
     11  157 1202 (S15E69 L= 20)   646     48        1 ( 138/1176- 178/1229) -20.6 E73-E65 S16-S13
     12  755 1189 (S03E19 L= 70)   268     29        0 ( 731/1165- 773/1219) -71.6 E20-E17 S04-S01
     13  744 1052 (S11E16 L= 73)   344    102        0 ( 702/1005- 777/1093)  71.6 E17-E14 S14-S08
     14 1401  891 (S02W25 L=114)  1121    236        4 (1364/ 775-1450/ 991) -36.0 W20-W31 S06-N03
     15  343  927 (S28E43 L= 46)   358    122        0 ( 310/ 866- 362/ 968) -63.4 E46-E40 S31-S25
     16 1725  736 (S00W53 L=142)   645    177        0 (1712/ 644-1736/ 816) -56.3 W48-W57 S05-N04
  

    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