Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  965 1788 (N52E23 L=222)   510    110        2 ( 914/1770-1015/1809)  26.6 E27-E17 N49-N53
      2  366 1457 (N18E55 L=190)  1276    420        5 ( 300/1334- 471/1550) -26.6 E64-E41 N11-N22
      3  809 1476 (N25E21 L=223)   218    166        0 ( 785/1401- 830/1542) -77.5 E23-E20 N20-N29
      4 1197 1354 (N23W06 L=251)   232     87        2 (1174/1305-1220/1394) -59.0 W04-W08 N20-N25
      5  209 1212 (N01E64 L=180)   453     80        0 ( 190/1186- 226/1237)  -8.1 E67-E60 N00-N01
      6 1553 1176 (N17W32 L=277)   246     87        1 (1524/1134-1575/1216) -68.2 W30-W34 N13-N18
      7 1341 1178 (N14W18 L=262)   227     88        1 (1301/1160-1370/1208)  45.0 W15-W19 N12-N16
      8 1385  944 (N00W23 L=267)   954    242        2 (1341/ 857-1435/1058) -74.7 W20-W27 S04-N07
      9  338  855 (S19E45 L=199)   419    132        0 ( 301/ 816- 393/ 901) -20.6 E48-E40 S21-S16
     10 1466  678 (S15W34 L=278)  1033    298        7 (1394/ 553-1523/ 776)  78.7 W31-W35 S23-S08
     11 1002  722 (S18W03 L=247)   250    102        0 ( 974/ 676-1042/ 749)  68.2 W01-W04 S21-S16
     12 1217  478 (S31W22 L=267)   324     58        1 (1185/ 453-1245/ 495)   0.0 W19-W24 S33-S30
  

    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