Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  533 1616 (N32E49 L=285)  1290    173        4 ( 500/1535- 586/1677) -59.0 E54-E45 N26-N36
      2 1531 1397 (N30W32 L=  6)   530    139        0 (1492/1330-1572/1452) -54.5 W29-W35 N26-N33
      3  643 1274 (N11E27 L=307)   665    200        2 ( 623/1189- 669/1366) -70.0 E30-E24 N06-N17
      4 1043 1209 (N12E01 L=333)   728    273        3 ( 990/1088-1078/1316) -61.7 E05-W02 N05-N18
      5 1348 1130 (N11W19 L=353)   627    184        2 (1304/1051-1389/1211)  70.0 W16-W21 N05-N16
      6  229 1004 (S10E58 L=276)  2600    368        5 ( 197/ 839- 265/1171)  81.5 E63-E53 S20-N00
      7 1877 1069 (N13W66 L= 40)   869     86        0 (1833/1045-1916/1093)  -4.1 W59-W73 N11-N14
      8 1446  548 (S24W36 L= 10)  2661    577        5 (1215/ 465-1678/ 627)  25.1 W21-W53 S31-S16
      9  946  519 (S32W01 L=335)   679    163        3 ( 869/ 505-1027/ 540)   5.2 E04-W07 S34-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