Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1250 1817 (N45W48 L=320)  2109    282        3 (1142/1748-1346/1893)   4.6 W32-W69 N41-N47
      2 1520 1423 (N14W40 L=313)  1356    266        2 (1420/1361-1634/1492) -50.2 W35-W46 N07-N19
      3  317 1472 (N47E48 L=224)   414     54        0 ( 287/1457- 352/1487) -11.3 E53-E43 N45-N48
      4  841 1429 (N31E01 L=271)   375     80        4 ( 814/1387- 868/1469)  59.0 E03-W00 N28-N33
      5 1056 1399 (N24W11 L=283)   244     63        2 (1020/1366-1092/1426) -11.3 W08-W13 N21-N25
      6  555 1320 (N31E23 L=249)  1395    281        2 ( 419/1296- 693/1354) -18.4 E34-E13 N27-N34
      7 1026 1251 (N16W05 L=277)   504    236        0 ( 933/1199-1129/1303)   4.4 E01-W12 N15-N17
      8  778 1017 (N08E14 L=258)   951    365        3 ( 639/ 980- 961/1058) -11.9 E22-E03 N04-N11
      9 1468  978 (S09W24 L=297)   554    119        1 (1427/ 917-1501/1036) -90.0 W22-W25 S13-S05
     10  441  999 (N14E36 L=236)   243     48        1 ( 412/ 981- 470/1015)   0.0 E38-E33 N13-N15
     11 1677  964 (S15W39 L=311)   275     42        0 (1656/ 936-1692/ 985) -90.0 W37-W40 S17-S13
     12  270  710 (N01E58 L=214)   726     83        3 ( 243/ 666- 305/ 748)  23.2 E61-E53 S01-N03
     13  665  651 (S10E30 L=243)   561    101        1 ( 624/ 601- 700/ 702) -80.5 E32-E28 S12-S06
     14  611  640 (S09E34 L=239)   243     66        0 ( 587/ 615- 630/ 666) -71.6 E35-E32 S10-S07
     15 1808  555 (S47W57 L=330)   598     45        0 (1779/ 534-1832/ 571) -14.9 W51-W66 S50-S44
     16 1361  535 (S33W09 L=282)   428     46        1 (1330/ 513-1396/ 559)  11.3 W07-W12 S34-S31
     17  387  420 (S18E63 L=209)   520     52        2 ( 365/ 395- 408/ 445)  71.6 E65-E59 S19-S16
     18  786  357 (S31E34 L=238)  1257    202        1 ( 687/ 301- 850/ 432)  13.2 E43-E26 S34-S26
     19  779  171 (S45E52 L=220)  1414    126        2 ( 720/ 132- 821/ 221)   6.0 E64-E43 S47-S40
  

    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