Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1797 (N30E53 L= 32)  1580    258        3 ( 603/1733- 662/1868) -54.0 E59-E48 N25-N36
      2 1623 1431 (N35W27 L=112)  2856    404        6 (1467/1222-1733/1590) -33.7 W13-W38 N24-N42
      3  287 1526 (N06E68 L= 17)  1089    236        2 ( 253/1489- 344/1566) -16.4 E77-E60 N03-N08
      4  433 1537 (N08E53 L= 32)  1262    201        0 ( 355/1509- 516/1563)  17.1 E59-E45 N05-N10
      5  423 1306 (S05E44 L= 41)  2016    359        2 ( 309/1129- 531/1452)  72.3 E48-E40 S17-N05
      6  940 1427 (N13E16 L= 69)   287     71        1 ( 895/1411- 985/1448)  45.0 E17-E13 N11-N15
      7 1316 1242 (N13W10 L= 96)  1902    321        3 (1164/1043-1417/1387) -30.6 E01-W21 N04-N17
      8  155 1034 (S26E62 L= 23)  1068    295        0 ( 103/1009- 216/1076)   0.0 E69-E55 S27-S23
      9 1575  655 (S11W44 L=129)   306     79        0 (1563/ 608-1592/ 695) -59.0 W42-W46 S13-S08
     10 1563  441 (S22W55 L=140)   723    140        0 (1504/ 425-1616/ 457)  24.0 W49-W59 S23-S19
  

    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