Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  995 1546 (N31W04 L=109)   287     56        2 ( 974/1515-1019/1577) -76.0 W02-W05 N29-N33
      2  262 1338 (N26E52 L= 53)  2123    381        2 ( 207/1194- 303/1447) -70.6 E60-E47 N16-N33
      3 1329 1250 (N09W21 L=126)   415    113        1 (1278/1196-1370/1291) -60.3 W17-W22 N05-N12
      4  207 1201 (N18E57 L= 49)   448     78        2 ( 193/1163- 225/1244)  68.2 E58-E55 N15-N20
      5 1779  748 (S25W52 L=157)   398    146        2 (1736/ 734-1821/ 767)  -7.1 W47-W55 S25-S23
      6  476  726 (S14E39 L= 67)   241     51        0 ( 452/ 711- 511/ 747) -33.7 E40-E36 S14-S12
      7 1635  533 (S38W42 L=148)   542     89        1 (1591/ 492-1681/ 569)  33.7 W38-W45 S40-S35
      8  405  479 (S28E52 L= 53)  1767    290        2 ( 349/ 412- 456/ 559)  36.9 E61-E48 S31-S22
  

    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