Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  375 1497 (N32E46 L=327)   637    135        0 ( 348/1424- 417/1548)  69.4 E47-E42 N27-N35
      2 1877 1335 (N13W68 L= 80)   782    138        1 (1837/1291-1902/1375) -39.8 W63-W71 N10-N15
      3  320 1251 (N16E47 L=325)   429     71        1 ( 299/1217- 338/1286) -53.1 E48-E45 N14-N18
      4 1074 1235 (N10W04 L= 17)   213     45        1 (1055/1212-1089/1270)  76.0 W02-W05 N08-N12
      5 1469 1240 (N08W29 L= 41)   276    100        0 (1423/1209-1512/1258) -31.0 W26-W31 N06-N09
      6  528 1188 (N11E30 L=342)   336     84        1 ( 488/1163- 573/1206)   0.0 E33-E27 N09-N11
      7  644 1113 (N06E23 L=350)   548    136        0 ( 588/1044- 676/1164) -58.0 E26-E20 N01-N09
      8  283 1044 (N04E50 L=322)   659    127        1 ( 260/ 985- 317/1084)  80.5 E52-E46 N00-N06
      9  383 1025 (N02E42 L=331)   411     99        1 ( 349/ 998- 427/1052) -18.4 E44-E38 N00-N03
     10 1819  933 (S11W55 L= 68)   310     51        0 (1793/ 913-1841/ 951)  11.3 W52-W57 S12-S10
     11 1660  684 (S26W43 L= 55)   535    138        0 (1617/ 621-1714/ 742)  49.4 W40-W46 S29-S22
     12 1190  447 (S40W08 L= 21)   264     48        0 (1159/ 433-1218/ 463)  14.0 W06-W10 S40-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