Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1187 1492 (N26W16 L= 50)   941    229        2 (1070/1440-1277/1522) -30.3 W08-W21 N21-N28
      2 1597 1380 (N15W42 L= 76)   322     79        2 (1566/1350-1627/1410) -53.1 W40-W43 N12-N17
      3  509 1344 (N23E31 L=  3)   977    222        3 ( 417/1281- 608/1392) -31.0 E39-E23 N18-N27
      4 1067 1376 (N20W06 L= 40)   200    108        0 (1010/1357-1109/1389) -18.4 W02-W08 N18-N20
      5 1348 1367 (N16W24 L= 58)   223     86        0 (1314/1346-1388/1385)   0.0 W21-W26 N15-N17
      6 1034 1156 (N07W02 L= 36)   178     94        1 ( 996/1121-1065/1198) -59.0 W00-W03 N04-N09
      7 1256 1097 (N01W14 L= 48)   224    146        1 (1198/1058-1302/1152)  23.2 W10-W17 N00-N03
      8 1894  981 (S11W65 L= 99)  1520    254        5 (1854/ 874-1940/1095) -52.4 W59-W73 S17-S04
      9 1700  934 (S13W44 L= 78)   282     61        1 (1669/ 910-1729/ 952) -11.3 W41-W46 S13-S11
     10 1500  710 (S25W29 L= 63)   502    215        2 (1474/ 624-1521/ 807)  90.0 W26-W30 S30-S18
     11  635  522 (S28E33 L=  1)  1333    260        1 ( 570/ 417- 703/ 629)  39.1 E42-E25 S34-S21
  

    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