Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1423 1836 (N40W69 L=283)   773     82        1 (1385/1811-1461/1851)  -7.1 W62-W79 N38-N41
      2 1401 1677 (N31W47 L=262)   778    196        1 (1314/1651-1493/1700) -26.6 W40-W53 N27-N34
      3 1071 1519 (N31W16 L=231)   615    228        0 ( 964/1495-1166/1544) -17.1 W09-W22 N27-N33
      4  745 1479 (N37E07 L=208)   974    295        1 ( 613/1434- 884/1526)  -5.4 E17-W04 N35-N38
      5 1450 1346 (N11W33 L=248)   757    210        0 (1395/1251-1496/1429) -90.0 W32-W35 N04-N17
      6 1886 1158 (S11W65 L=280)  1238    140        2 (1847/1117-1948/1213) -21.8 W61-W77 S13-S06
      7  364  874 (N09E45 L=169)  1745    522        0 ( 211/ 695- 535/ 997) -46.4 E58-E37 S03-N18
      8 1733  637 (S39W43 L=258)  2705    376        3 (1639/ 516-1834/ 799)  35.0 W33-W54 S45-S30
      9 1360  734 (S21W13 L=227)   572    136        0 (1316/ 684-1387/ 797)  39.8 W09-W15 S23-S17
     10  665  421 (S23E40 L=174)   532    190        0 ( 644/ 338- 691/ 507)  45.0 E45-E34 S29-S18
     11  507  407 (S21E53 L=161)  1072    193        2 ( 477/ 321- 541/ 497)  45.0 E58-E46 S27-S15
     12 1220  278 (S50E10 L=205)  1080    155        3 (1147/ 236-1300/ 325)   0.0 E16-E01 S52-S47
  

    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