Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1629 (N42W22 L=113)   271     60        2 (1048/1599-1083/1659)  45.0 W19-W24 N40-N44
      2  446 1434 (N45E36 L= 55)   391     90        1 ( 403/1413- 491/1459)   0.0 E40-E30 N44-N46
      3  433 1316 (N37E36 L= 55)   315     60        1 ( 402/1295- 467/1330)   0.0 E39-E32 N36-N38
      4  251  924 (N16E57 L= 34)  1567    256        4 ( 194/ 827- 305/1048) -61.9 E62-E52 N09-N24
      5  186  914 (N16E65 L= 26)  1088    164        3 ( 170/ 830- 207/ 992)  84.3 E68-E63 N10-N20
      6 1618  830 (S19W33 L=124)   418    180        0 (1602/ 742-1645/ 904)  56.3 W30-W36 S24-S15
      7 1366  701 (S20W13 L=104)  1584    527        1 (1186/ 533-1511/ 880)  25.7 E02-W25 S26-S12
      8 1091  673 (S15E04 L= 87)   755    280        2 (1031/ 582-1156/ 793)  45.0 E09-W02 S19-S08
      9  256  612 (S03E68 L= 23)   681     65        3 ( 234/ 576- 283/ 649) -53.1 E71-E64 S05-N00
     10 1061  476 (S27E12 L= 79)   803    170        1 (1014/ 439-1125/ 547) -77.5 E14-E07 S30-S21
     11 1399  495 (S35W11 L=102)   393     91        0 (1360/ 461-1436/ 535)  23.2 W07-W14 S36-S33
     12  720  252 (S36E48 L= 43)  2137    230        2 ( 667/ 157- 765/ 352)  36.9 E68-E42 S44-S26
  

    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