Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1558 1701 (N32W60 L= 25)   846    162        2 (1520/1678-1605/1727) -33.7 W56-W63 N29-N33
      2 1020 1587 (N41W15 L=340)   449     86        1 ( 992/1547-1041/1625)  68.2 W12-W17 N38-N43
      3 1148 1576 (N37W23 L=349)   245     58        1 (1112/1561-1182/1589)  -9.5 W20-W26 N36-N38
      4 1765 1419 (N11W63 L= 28)   920    219        0 (1744/1368-1806/1468)  26.6 W59-W71 N06-N12
      5 1015 1423 (N31W09 L=334)   347    106        2 ( 961/1403-1069/1446) -29.7 W05-W12 N28-N33
      6 1683 1367 (N10W52 L= 17)   453    143        1 (1629/1324-1730/1396) -24.0 W48-W57 N07-N13
      7 1299 1244 (N14W21 L=346)   245     46        1 (1270/1216-1322/1267) -76.0 W19-W22 N11-N15
      8 1004 1185 (N17W02 L=327)   204    121        0 ( 942/1166-1066/1196) -33.7 E01-W05 N14-N18
      9 1323 1054 (N02W18 L=343)   398     93        0 (1297/1005-1359/1102) -80.5 W16-W19 S01-N05
     10  275  849 (N09E55 L=270)   421    105        1 ( 253/ 815- 310/ 904) -59.0 E56-E52 N07-N12
     11 1428  710 (S21W18 L=343)   433    141        0 (1380/ 640-1478/ 770)  29.1 W13-W22 S23-S17
     12 1234  422 (S35E00 L=325)   515    105        0 (1195/ 380-1281/ 475)  35.5 E04-W04 S38-S32
  

    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