Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1795 1430 (N45W49 L=138)   844    110        1 (1747/1394-1836/1453)  14.0 W41-W57 N41-N47
      2  584 1204 (S08E30 L= 59)   334     63        1 ( 557/1174- 615/1239)  68.2 E31-E28 S10-S05
      3 1760 1064 (N18W45 L=135)   596    151        1 (1720/ 997-1811/1118) -16.7 W41-W51 N16-N21
      4  609  986 (S20E23 L= 66)   447    133        0 ( 553/ 940- 660/1021)  49.4 E26-E20 S23-S16
      5  433  986 (S24E36 L= 53)   536    192        1 ( 357/ 951- 513/1007)  13.0 E42-E29 S26-S21
      6  327  775 (S40E44 L= 45)   799    135        0 ( 281/ 728- 377/ 843)  45.0 E50-E38 S43-S33
      7 1327  517 (S25W36 L=125)  1311    214        4 (1255/ 459-1433/ 559)  63.4 W29-W40 S29-S19
      8 1033  324 (S47W30 L=119)  1183    203        1 ( 997/ 271-1106/ 413) -15.5 W21-W39 S49-S41
      9 1195  373 (S37W37 L=126)   456     56        1 (1161/ 345-1220/ 404) -15.9 W32-W39 S39-S35
  

    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