Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1483 1620 (N31W42 L=202)   299     99        0 (1462/1582-1499/1653)  45.0 W38-W44 N28-N33
      2 1698 1396 (N14W51 L=212)   630    179        1 (1659/1332-1745/1434)  56.3 W47-W56 N10-N16
      3  160 1144 (N17E62 L= 98)  4094    585       14 (  77/1026- 225/1258)  31.0 E77-E54 N10-N23
      4  369 1223 (N19E41 L=119)   305     60        1 ( 341/1200- 397/1241) -26.6 E43-E39 N17-N20
      5  448 1087 (N10E36 L=124)   381    103        0 ( 405/1039- 489/1131) -50.2 E38-E33 N06-N12
      6  616  983 (N02E25 L=135)   478    106        1 ( 585/ 920- 640/1034)  51.3 E28-E23 S01-N04
      7 1114  808 (S14W03 L=163)   282    113        0 (1061/ 788-1181/ 841)   8.1 W00-W07 S15-S13
      8 1076  651 (S24E01 L=159)   520    215        0 ( 984/ 619-1190/ 708)   8.7 E07-W06 S25-S21
      9  993  426 (S38E11 L=150)   311     62        0 ( 968/ 398-1015/ 459)  38.7 E13-E08 S39-S35
     10  599  397 (S35E41 L=119)   433    117        0 ( 582/ 354- 618/ 444)  50.2 E44-E39 S37-S31
     11  912  370 (S41E19 L=142)   335     54        1 ( 888/ 350- 936/ 401)  45.0 E20-E16 S42-S38
  

    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