Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  846 1811 (N44E41 L= 53)  2955    313        2 ( 793/1696- 931/1942) -22.2 E75-E25 N36-N56
      2  340 1458 (N09E58 L= 35)  3811    366        8 ( 241/1247- 428/1657) -77.5 E68-E55 S05-N22
      3  567 1447 (N12E39 L= 54)   298     63        2 ( 545/1425- 597/1473) -18.4 E41-E35 N11-N13
      4  509 1267 (N01E37 L= 56)   241     65        0 ( 489/1229- 522/1298) -45.0 E38-E35 S01-N02
      5 1824 1222 (N28W55 L=148)   756    146        1 (1790/1163-1861/1270) -45.0 W50-W59 N25-N31
      6  159 1179 (S11E69 L= 24)  1477     91        4 ( 118/1145- 192/1209)  -8.1 E78-E63 S12-S08
      7  996 1062 (S01E02 L= 91)   373     78        1 ( 974/1028-1019/1102)  90.0 E03-E01 S03-N01
      8 1419  913 (S00W26 L=119)   683    135        1 (1375/ 837-1448/ 967) -29.7 W21-W28 S04-N02
      9 1480  589 (S17W40 L=133)   333      9        0 (1460/ 567-1503/ 610) -14.0 W37-W41 S18-S15
     10 1758  496 (S15W74 L=167)   693     79        0 (1721/ 474-1797/ 519)  29.1 W69-W79 S17-S12
     11  990  247 (S54W26 L=119)   317     67        0 ( 979/ 214-1004/ 283) -40.6 W22-W29 S56-S50
  

    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