Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1192 1495 (N26W16 L= 49)   972    216        0 (1077/1443-1285/1523) -24.8 W08-W21 N21-N28
      2  511 1348 (N24E31 L=  2)   959    212        2 ( 423/1291- 610/1393) -28.1 E38-E23 N19-N27
      3 1073 1379 (N20W07 L= 40)   239    103        0 (1016/1360-1119/1390) -18.4 W03-W09 N18-N20
      4 1132 1240 (N11W09 L= 42)   271    107        1 (1110/1210-1174/1271) -53.1 W07-W10 N08-N12
      5 1042 1158 (N07W02 L= 36)   216    112        0 (1001/1120-1077/1207) -51.3 W00-W04 N04-N09
      6  806 1081 (N04E12 L= 21)   215     90        0 ( 796/1029- 819/1123)  78.7 E13-E11 N01-N06
      7 1868 1017 (S09W61 L= 95)   887    172        0 (1830/ 943-1892/1096)  56.3 W56-W64 S13-S04
      8 1706  935 (S12W45 L= 78)   253     53        2 (1680/ 911-1735/ 954)   0.0 W42-W47 S13-S11
      9 1387  915 (S11W21 L= 54)   214     48        1 (1365/ 883-1409/ 935) -71.6 W19-W22 S12-S09
     10 1512  674 (S27W30 L= 63)   370     93        1 (1499/ 623-1528/ 725) -74.1 W28-W31 S30-S23
     11  209  680 (S13E67 L=327)   492     91        2 ( 182/ 648- 232/ 719) -78.7 E69-E64 S15-S10
     12  641  523 (S28E33 L=  1)  1475    259        2 ( 576/ 416- 710/ 631)  37.4 E42-E24 S34-S21
  

    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